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Abstract. - We use a many-body, atomistic empirical pseudopotential approach to predict the 
multi-exciton emission spectrum of a lens shaped InAs/GaAs self-assembled quantum dot. We 
discuss the effects of (i) The direct Coulomb energies, including the differences of electron and 
hole wavefunctions, (ii) the exchange Coulomb energies and (iii) correlation energies given by a 
configuration interaction calculation. Emission from the groundstate of the N exciton system 
to the N — 1 exciton system involving eo — > ho and ei — > hi recombinations are discussed. A 
comparison with a simpler single-band, effective mass approach is presented. 

High-resolution single-dot spectroscopy (fl-pl of InAs/GaAs self-assembled quantum dots 
shows that as the excitation intensity is increased, thus loading more excitons into the dots, 
new emission lines appear both to the red and to the blue of the fundamental emission line 
observed at low excitation power. "State filling" effects, leading to the recombination of high 
energy electron-hole pairs, cannot explain the red-shifted emission lines, nor the fact that the 
number of lines exceeds the number of allowed single-particle transitions. In this letter we 
present a theory of self-assembled semiconductor quantum-dots, based on a pseudopotential 
many-body expansion that demonstrates that it is multi-exciton transitions that are responsi- 
ble for this complex observed spectral structure. We isolate and clarify three distinct physical 
effects; (i) electron-hole wavefunction asymmetry, leading to a blue shift of the fundamen- 
tal exciton transition as the number of spectator excitons loaded into the dot increases, (ii) 
electron-electron and hole-hole exchange interactions which red shift all even multiexciton de- 
cays (biexciton, four-exciton) and split the odd multiexciton decays (tri-exciton, five-exciton) 
into multiple-lines and (iii) correlation effects which red shift the biexciton leading to its 
binding with respect to the monoexciton. 

We will first describe the qualitative picture of multi-excitons and then describe a quan- 
titative model. The essential physics of such transitions can be understood by considering 
what happens to the ground-state recombination of the lowest electron level, eo, and the low- 
est hole level, ho, if other electrons and holes are present in the dot as "spectators". The 
schematic figures in the center of Fig. |l| depict the fundamental eo — ho recombination in the 
presence of to 5 "spectator" electron-hole pairs (we assume here that all levels are spatially 
non-degenerate Q). We distinguish here two exciton series; (i) when the initial number of 
excitons, N is even, the initial electron configuration is "closed shell", e.g. (eg)(/io) f° r N = 2, 
whereas (ii) when N is odd, the initial configuration has an open shell both in the electron 
and in the hole manifold, e.g. (efyeKhfyhl for N = 3, (parentheses are used for the closed 
shell orbitals). The distinction between the "closed shell" and "open shell" multiexciton is 
important, since in the initial state of the A^=even, "closed shell" series the spectator levels 
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can be occupied in only a single unique manner by the spectator spins, while in the ./V=odd, 
"open-shell" series, many spin arrangements are possible in the initial and final states. This 
will lead to a large number of exciton lines. 

To understand qualitatively the effects of spectator electrons and holes, let us consider 
the Hartree-Fock (HF) energy of a single configuration of electrons and holes. We denote 
the direct Coulomb interaction between carriers in levels i and j by Jij and the exchange 
interactions as K%j . The recombination energy of the mono-exciton is given by 

^e^hg = ( £ e ~ e h ) ~ JeaM ' (1) 

where e eo and e^ are the single particle levels and -E^o denotes the energy associated with 
decaying from a single exciton to the groundstate. Neglecting the electron-hole exchange 
(NB. it is included in the calculation) the eo — ho recombination energy in the presence of N s 
electrons and holes is 

(2) 

where e s and h s are spectator electrons and holes, such that Yle* ~ = 1. We see from 
Eq.(^) that the eg — ho recombination energy is shifted with respect to the fundamental exciton. 
This shift has two sources indicated by the two bracketed terms in Eq.(||). First, the Coulomb 
shift, which vanishes if the electrons and holes, eo and ho have the same wavefunctions, i.e. if 
Je e s = Je s h aim Jh h s = Je h s - This "Coulomb shift" , S^ ! °_^ l N _ 1 thus reflects the difference 
in the electron and hole wavefunctions. It vanishes artificially in single-band effective-mass 
calculations that use an infinite well depth, such as those in Refs. Second, there is 

an exchange shift, S'^^ L N _ 1 , given by the second term in brackets in Eq.fQ). This exchange 
shift is familiar from theories of band gap renormalization 0] where the existence of high 
carrier densities during high power photoexcitation act to reduce the band gap. In addition, 
since the exchange interaction depends on the spin orientation of the carriers, the exchange 
contribution [second term in Eq.(H)] can split the excitonic transitions. 

In the iV=even, "closed shell" series, the initial state contains no open shells while the final 
states contains one open shell. This results in one- and four-fold degeneracies for the initial 
and final states, and hence 4 transitions. These 4 transitions are split by the small electron- 
hole exchange interaction. In contrast, for the iV=odd transitions, the initial and final states 
contain one and two open shells respectively, resulting in four- and 16-fold degeneracies and 
hence a total 64 possible transitions. Different alignments of the spins produce splittings of 
the different transitions resulting from electron-electron and hole-hole exchange interactions. 
In summary, within the HF approximation the presence of "spectator excitons" will shift the 
e o — ho transition in the 7V=even series and shift and split the eo — ho transition in the ./V=odd 
series. 

The above treatment neglects the effects of correlation. These can be conveniently intro- 
duced by considering configuration-interaction (CI) effects [Q. In other words, instead of 
evaluating the energy [Eq. ([[])] of the e - ho monoexciton by considering only the 4 spin ar- 
rangements consistent with a single orbital configuration, ej/ij, we allow the presence of other 
orbital configurations such as ej/ij, t\h\ and &\h\. This configuration interaction approach 
can shift the Hartree-Fock transitions of Eqs.(|l|) by a "correlation shift" , 6^ \ , which includes 
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both self-consistent adjustments of the single particle orbitals and correlation. This will, for 
example, provide additional binding to the biexciton with respect to the monoexciton 

p2— >1 _ pl^O [ jee , jhh _ n jeh 1 _ ?CI (o\ 

^eoho ^eoho ~ l J e ,e T J h ,h z<J e ,/i J u e h • \° ) 

Correlation effects can also add new transitions to the spectrum, due to mixing of new config- 
urations into the Hartree Fock configuration. Thus, CI will, in general, also alter the intensity 
of the HF transitions. 

The purpose of this letter is to make a realistic prediction of the "Coulomb shift" , S Coul , 
the "exchange-splitting', S Exch , ' and "correlation shift/splitting", 5 CI , due to the presence 
of spectator excitons in a self-assembled semiconductor quantum dot. We adopt the exper- 
imentally determined lens shaped dot with a base of 250 A and a height of 35 Aand a PL 
peak at 1.1 eV. The optical and electronic properties of this dot geometry have been exten- 
sively studied, see Refs. jll],[l^] and references therein. Recently, alternative dot geometries 
for uncapped dots have been proposed, based on {136} facets fl3|| . However, there is strong 
evidence of Ga in-diffusion during the capping process which acts to produce a lens shaped 
geometry. The single-particle bound states are calculated using an empirical pseudopoten- 
tial Hamiltonian || . Recent developments in these pseudopotentials enable the full inclusion 
of the effects of strain, multi-band couplings, band non-parabolicity and spin-orbit coupling 
in the single particle Schrodinger equation. In Ref. H we demonstrate that by accurately 
fitting the bulk band structure, effective masses and deformation potentials we are able to 
use these pseudopotentials to obtain excellent agreement with a wide range of optical and 
electronic properties for the lens shaped dots discussed in this paper. These comparisons were 
not possible earlier when the shape of the dot and its composition profile were unknown. In 
Ref. H , we model the shape and composition profile of the dot and then compare calculated 
and measured single exciton energies, inter-band energy spacings, electron and hole binding 
energies and wetting layer energies, finding excellent agreement. For completeness, we treat 
here the same lens-shape, alloyed dot with a 1.1 eV PL peak. Unfortunately, although nu- 
merous experiments were conducted on this dot, no multi-exciton spectra were taken since 
the lowest PL is outside the range of conventional CCD detection equipment, so quantitative 
comparison awaits a future measurement. 

Having obtained the single-particle levels (eo,ei..., ho, hi...) from the pseudopotential 
method, we then calculate numerically the screened Coulomb and exchange integrals 

J J e(ri - r 2 )|ri - r 2 | 

_ r I iSMfM^MliM dridr2 , (4) 

J J e(ri - r 2 )|ri - r 2 | 

where e is a size dependent, phenomenological, screened dielectric function |L4|. Our exchange 
automatically includes both short and long range components |l5| ]. In the configuration- 
interaction approach, we expand the wavcfunction of the A^-exciton, <]/, as a linear combination 
of Slater determinants, $, obtained by exciting N electrons from the valence band to the 
conduction band. For example, in the N = 2 (biexciton) case we have: 

*(ri, • • ■ ,r M ) = y2 A w',cc'®vv',cc>(ri,--- ,*m) , (5) 



Here M is the total number of electrons in the quantum dot, v, v' denote the unoccupied 
valence band single particle states (holes) and cc' denote the conduction band states occupied 
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by the N excited electrons. To obtain our many-body states, we diagonalize the many-body 
Hamiltonian in the basis of $. 

To clarify and isolate the physical factors contributing to multi-exciton effects we solve the 
problem in a series of steps. First, we neglect configuration interaction effects, treating only 
single-configurations and set all the exchange integrals to zero, (5 exch — 0). In this (e + J) 
approximation we see the effect of the Coulomb "chemical shift" , 8 : . Second, we will still 
use a single-configuration, but include exchange integrals. In this (e+J+K) approximation we 
will see the added effects of carrier-carrier exchange, 8 exch . Finally, the effects of correlation 
are included in the CI calculation. 

Figure ^(a) shows our calculated energies associated with the fundamental recombination 
of an electron and hole, e — > h a in the presence of to 5 spectator excitons. These are denoted 
as the 1 — > to 6 — > 5 transitions. The N excitons that form the initial state of each transition 
are assumed to occupy the groundstate configuration as predicted by the Aufbau principle. 
A schematic energy level scheme in the center of each figure, shows this initial occupation, 
where the vertical solid line marks the recombination taking place. The spectrum shown in 
each panel is obtained from a sum of Gaussians (0.1 meV width) whose means represent the 
transition energies and heights are proportional to the calculated dipole transition element. 
The zero of energy is taken as the fundamental exciton, e eo — e/j — Je h - The multiplicity of 
each individual group of transitions in the (e + J) and (e + J + K) approximations are marked 
in black and red. 

The effect of direct Coulomb interactions: The red lines in Fig. |l](a) show the transition 
energies in the (e + J) approximation which include only single particle and direct Coulomb 
energies. Inspection of these lines shows that: 

(i) All the observed Coulomb shifts are relatively small, (6 \ Coul ~ 2 meV) and result in a 
blue shift of the transition with respect to the fundamental transition. For example, at this 
level of approximation, the biexciton is "unbound" with respect to two single excitons, i,e. 

(ii) As all the transitions involve the same single particle eo — * ho recombination, they all 
have the same oscillator strength. 

(iii) All the transitions show only a single degenerate line as neither the initial or the final 
state exhibit any exchange splittings in this approximation. 

The effect of exchange interactions: The black lines in Fig. |l|(a) show the transition 
energies in the (e + J + K) approximation which includes single particle, direct and exchange 
Coulomb energies. Inspection of these lines shows that: 

(i) The 1 — v and 2 — ► 1 transitions contain only one or zero spectator excitons and hence 
no electron-electron or hole-hole exchange takes place. Therefore, only the small electron-hole 
exchange fine splitting is observed. We therefore classify these transitions as essentially 4-fold 
multiplets. Due to the lack of correlation the biexciton is still unbound with respect to two 
single excitons as in the (e + J) approximation. 

(ii) For even — > odd there are 4 possible transitions. In each case, 2 — ► 1, 4 — ► 3 and 
6 — > 5 we observe only one 4- fold multiplet. The fine splitting within this group of tran- 
sitions arises from electron-hole exchange splittings of the 4 final states. The initial states 
contain no unpaired electrons and holes and hence exhibit no exchange splittings. As the 
number of spectator excitons increases, there is a red shift of these transitions due to ex- 
change interactions. The exchange energy shifts from Eq.(|J) are (Jffjf = [K eoei + Kh^] and 
fiexch _ [K eoe2 + K eoei + K hoh2 + K hohl ]. 

(iii) For odd — > even with N > 3 there are 64 possible transitions. In both cases: 3^2 
and 5 — > 4 we see 6 groups of transitions with multiplicities of 4:4:8:12:12:24. The splitting 
between the 6 groups arises from electron-electron and hole-hole exchange splittings between 
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the final states which contain two unpaired electrons and holes. The 6 groups of transitions in 
3^2 and 5-^4 span 23 and 22 meV. This energy span reflects the span of the 16 eigenvalues 
of the 16x16 matrix generated by the different spin occupations of the final states. 

(iv) The exchange interaction alters the oscillator strength of the transitions so that they 
are not all identical as in the (e + J) approximation. 

Comparison with previous calculations: Previous calculations of multi-excitons in quantum 
dots include the works of Hu |T(|],Takagahara JL7[ , Barenco and Dupertuis 0, Dekel et. al. 
[jl],^], Landinet. al. || and Hawrylak Refs. |L|,p|,|7| JT^ , p7[ and p| adopt single-band effective- 
mass models with either an infinite potential barrier pj, |2|, [7|, pl| |l7| or a parabolic potential J8j, 
both of which artificially force the electron and hole wavefunctions to be identical. Therefore, 
in all these calculations the Coulomb shift, S Coul , is zero and the electron-hole exchange 
vanishes, K ei h- — 0. To obtain realistic single particle energy spacings in Refs. Jl|,|2[0] an 
unrealistic cuboidal shape had to be assumed and the size of the dots was treated as adjustable 
parameters. By choosing different lengths for all three dimensions both the measured s-p and 
p-p splittings can be reproduced. The two-dimensional parabolic potential adopted in Ref. || 
can also be adjusted to reproduce the correct s-p splitting, but will always produce degenerate 
p states. It has recently been shown that only models which can split this p level degeneracy 
can provide a realistic interpretation of experimental results. In addition to approximating 
the single particle states, Refs. Ij],|^,ffl,[|,[l6 , 17 neglect the effects of strain and the spin-orbit 
interaction. 

To calculate the exchange and correlation contribution to the excitonic energies Ref. jl6| 
uses path integral quantum Monte Carlo techniques which provide an exact (to within statis- 
tical error) determination of the exchange and correlation energy. However, quantum Monte 
Carlo methods are currently restricted to single band Hamiltonians and cannot therefore pre- 
dict the Coulomb shift, S Coul , discussed above. Ref. || adopts a limited basis CI to estimate 
correlation energies. Refs. PrPt^Hl^] use only a single configuration approach which does not 
include correlation effects. 

To assess the above approximations used in [[j} [?], ||, [l(| we show in Fig. |l](b) a repeat of 
our calculations for the transition energies within the (e + J) and (e + J + K) approximations 
applying the assumptions adopted in Refs. 0[|[7|,||, namely J eiej = J} uhj = J ei hj,K eiej = 
Kh.hj , K eihj = and A so = 0. We observe that: 

(i) Within the (e + J) approximation, the chemical shift, S Coul is zero by definition so 
all the red lines lie on the zero of energy, (ii) Within the (e + J + K) approximation for 
even — ► odd, there are 4 possible transitions that are exactly degenerate as = 0. (iii) 

For odd — > even, there are 64 possible transitions, split into 4:24:36 multiplets, compared to 
the 4:4:8:12:12:24 multiplets obtained in the pseudopotential calculations. The reduction in 
the number of multiplets arises form the assumptions ip Ci = iphi and A = 0. Other than 
these changes we find that many of the qualitative feature noted in the calculations of Dekel 
et. al. are retained in the pseudopotential description. 

The effect of CI interactions: In Fig. || we contrast the transition energies from our 
pseudopotential calculations within the (e + J + K) approximation (black lines) with those 
from a CI calculation (red lines) which includes single particle, direct and exchange Coulomb 
and "correlation" effects. 

In the CI calculations we expand the many-body wavefunction in a basis of Slater deter- 
minants constructed from all possible orbital and spin occupations of the lowest 10 (including 
spin degeneracy) electron and lowest 10 hole single particle levels. For example, the biexci- 
ton basis contains 10 C2. 10 C2 = 2025 Slater determinants, and the 3-exciton basis contains 
10 C3. 10 C3 = 14400 Slater determinants. This basis neglects the contributions for higher lying 
bound states and continuum states. To investigate the effects of adopting this limited CI ba- 
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sis |l8| we have compared the results of diffusion quantum Monte Carlo calculations (DMC) 
and CI calculations for the exciton correlation energy in a model, single band system with 
equivalent size, band offsets and number of bound states. We find that our CI calculations 
retreive approximately 50-60% of the 5-6 meV of correlation energy obtained in the DMC 
calculations and can therefore be used as a lower bound for the effects of correlation. Using 
these calculations we find that 

(i) Correlation effects lower the energy of both the initial and final state of a transition. 
Our calculated correlation varies from 2-3 meV for a single exciton to 10 meV for multiple 
excitons. For all the transitions shown here, the correlation shift, S CI , for the initial state 
with N excitons is greater than that for the final N — 1 exciton state. Therefore, all the CI 
transition peaks are red shifted with respect to those from the (e + J + K) approximation, (ii) 
This (S CI ) red shift is larger for the 2 — > 1 transition than for the 1 — > transition and is able 
to overcome the Coulomb blue shift and "bind" the biexciton. (iii) As the number of spectator 
excitons increases, the difference in the red shift for the initial and final states decreases, so 
that the red shift of the transition energy decreases. For 3^2 and 4^3 transitions the 
red shifts of the transitions rapidly decrease, (iv) The mixing of configurations within the CI 
results in both additional peaks in the CI spectra and changes in the relative magnitude of 
the peaks, e.g. for 4^3 the single (e + J + K) peak is split into two equally strong multiplets 
of peaks in the CI spectra. The additional peaks result from the mixing in of excited states 
within the CI calculation. For example, the additional CI peaks blue shifted from the main 
peaks in the 4—^3 spectra result from excited states configurations mixed into the N = 3 
exciton. 

In conclusion, we present results of the first pseudopotential, many-body calculation of 
multi-exciton states within InAs/GaAs quantum dots. We are able to isolate the effects 
of the direct and exchange Coulomb interactions and correlation on the energies of N — > 
N — 1 excitonic transitions. We find that direct Coulomb energies introduce small blue shifts. 
Electron-electron and hole-hole exchange splittings which are responsible for the majority of 
the observed structure, introduce both red shifts and splittings. Correlation effects red shift 
all transitions and change the relative energies of transitions (e.g. bind the biexciton). 
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Fig. 1 - Energy of eo — * ho, recombinations in the presence of to 5 spectator excitons. Energies 
calculated in the (e + J) (red) and (e + J + K) (black) approximations are shown, (a) Shows our 
pseudopotential calculations and (b) shows calculations with the assumptions from Refs. The 
multiplicities of each line is labelled. 
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Fig. 2 - Energy of eo — * ho, recombinations in the presence of to 3 spectator excitons. The black 
and red lines show energies calculated in the (e + J + K) and CI approximations. 



